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ABSTRACT 

The  secondary  photometric  standard  star  number  79  for  the  FS  Aur  field  (Henden  &  Honeycutt  1997),  designated 
as  [HH97]  FS  Aur-79  (GSC  1874—399),  is  a  short-period  (0.2508  days)  eclipsing  binary  whose  light  curve  is  a  com¬ 
bination  of  the  1 3  Lyr  and  BY  Dra  type  variables.  High  signal-to-noise  ratio  multicolor  photometry  was  obtained  using 
the  US  Naval  Observatory  1  m  telescope.  These  light  curves  show  asymmetry  at  quadrature  phases  (the  O’Connell 
effect),  which  can  be  modeled  with  the  presence  of  starspots.  A  low-resolution  spectrum  obtained  with  the  3.5  m 
Wisconsin-Indiana-Yale-NOAO  telescope  at  orbital  phase  0.76  is  consistent  with  a  spectral  type  of  dK7e  and  dM3e. 
A  radial  velocity  curve  for  the  primary  star  was  constructed  using  24  high-resolution  spectra  from  the  9.2  m  Hobby- 
Eberly  Telescope.  Spectra  show  Ha  and  H/3  in  emission  confirming  chromospheric  activity  and  possibly  the  presence 
of  circumstellar  material.  Binary  star  models  that  simultaneously  fit  the  U,  B,  V,  R ,  and  radial  velocity  curves  are  those 
with  a  primary  star  of  mass  0.59  ±  0.02  Mg,  temperature  4100  ±  25  K,  and  mean  radius  0.67  R®,  just  filling  its  Roche 
lobe,  and  a  secondary  star  of  mass  0.31  ±  0.09  M0,  temperature  3425  ±  25  K,  and  mean  radius  0.48  R0,  just  within  its 
Roche  lobe.  An  inclination  angle  of  83°  ±  2°  with  a  center-of-mass  separation  of  1.62  R  .  is  also  derived.  Starspots, 
expected  for  a  rotation  period  of  less  than  1  day,  had  to  be  included  in  the  modeling  to  fit  the  O’Connell  effect. 
Key  words:  binaries:  close  —  binaries:  eclipsing  —  binaries:  spectroscopic  —  stars:  late-type 


1.  INTRODUCTION 

Henden  &  Honeycutt  ( 1997)  designated  this  star  as  the  79th 
standard  star  in  the  field  of  cataclysmic  variable  star  FS  Aur;  this 
star  is  also  cataloged  GSC  1874—399  in  the  Guide  Star  Catalog 
and  misreferenced  as  [HH95]  FS-79  by  S1MBAD,1 2 3  although  it 
should  be  [HH97]  FS-79.  In  this  paper  we  refer  to  this  star  as  FS 
Aur-79.  The  discovery  that  this  star  was  actually  an  eclipsing 
binary  was  described  by  Robertson  et  al.  ( 2004).  They  found  that 
the  orbital  period  was  approximately  0.25  days  and  that  the  tem¬ 
peratures  of  the  two  stars  were  not  equal,  indicating  that  the  two 
stars  were  not  in  contact.  Their  low  signal-to-noise  ratio  (S/N) 
BVR  photometry  from  2003-2004  enabled  them  to  determine 
preliminary  system  parameters  such  as  mass  ratio,  inclination 
angle,  effective  temperatures,  and  filling  factors.  This  indicated 
a  near-contact  binary  system  with  late-type  components. 

Good,  accurate,  fundamental  parameters  for  late-type  main- 
sequence  stars  are  based  on  only  a  handful  of  double-lined  spec¬ 
troscopic,  totally  eclipsing  binaries  with  K  and  M  components. 
Therefore,  there  are  considerable  gaps  in  our  interpolations  and 
knowledge.  This  is  also  true  when  determining  effective  temper¬ 
ature  for  cooler  M  dwarfs,  as  the  complexity  of  molecules  present 
in  their  atmospheres  increases  with  later  spectral  types.  Near¬ 

1  Michelson  Postdoctoral  Fellow. 

2  Also  at  NVI,  Inc.,  Greenbelt,  MD  20770,  USA. 

3  The  SIMBAD  database  is  maintained  by  the  CDS,  Strasbourg,  France. 


contact  and  contact  systems  with  these  very  late  components  are 
extremely  rare  and  are  therefore  interesting  for  what  they  can  tell 
us  about  binary  stellar  formation  and  evolution. 

To  better  constrain  the  system  parameters  we  obtained  high 
S/N  UB  VR  photometry  and  spectra  to  determine  the  radial  veloc¬ 
ity  curves.  In  §  2  we  describe  how  the  photometry  and  spectroscopy 
were  obtained  and  reduced.  In  §  3  we  present  the  photometric 
and  spectroscopic  analysis.  Section  4  is  a  discussion  of  our  re¬ 
sults,  and  §  5  is  a  summary  of  our  conclusions. 

2.  OBSERVATIONS  AND  REDUCTIONS 
2.1.  Photometry 

2.1.1.  Center  for  Backyard  Astrophysics  Photometry 

After  the  initial  discovery  of  this  eclipsing  variable  star 
(Robertson  et  al.  2004),  additional  light  curves  were  collected  by 
T.  C.  on  15  nights  from  2003  November  through  2004  February. 
The  interest  of  the  Center  for  Backyard  Astrophysics  (CBA)  in 
the  cataclysmic  variable  FS  Aur,  very  near  our  newly  discovered 
variable,  also  helped  spark  a  professional-amateur  observational 
campaign  by  some  members  of  the  CBA  to  collect  differential 
photometry  for  obtaining  times  of  eclipse  minima.  Data  were 
gathered  by  T.  C.  (CBA-Arkansas),  T.  Krajci  ( 2006;  CBA-New 
Mexico),  R  de  Ponthiere  (CBA-Lesve),  and  J.  Foote  (CBA-Utah) 
on  45  nights  from  2005  October  to  2006  January.  Details  of  the 
observatories  and  their  telescope/detector  configurations  can 
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TABLE  1 

CBA  Times  of  Minimum 


HJD 

Phase  Cycle 

Filter 

Observatory 

2,452,963.74445 . 

0.0 

R 

CBA-Arkansas 

2,452,963.86918 . 

0.5 

R 

CBA-Arkansas 

2,452,964.74756 . 

4.0 

R 

CBA-Arkansas 

2,452,964.87062 . 

4.5 

R 

CBA-Arkansas 

2,452,965.75061 . 

8.0 

R 

CBA-Arkansas 

2,452,966.00206 . 

9.0 

R 

CBA-Arkansas 

2,452,968.76013 . 

20.0 

R 

CBA-Arkansas 

2,452,971.89528 . 

32.5 

R 

CBA-Arkansas 

2,452,972.77382 . 

36.0 

R 

CBA-Arkansas 

2,452,984.68797 . 

83.5 

V 

CBA-Arkansas 

2,452,984.81383 . 

84.0 

V 

CBA-Arkansas 

2,452,985.68966 . 

87.5 

V 

CBA-Arkansas 

2,452,985.81665 . 

88.0 

V 

CBA-Arkansas 

2,452,988.70157 . 

99.5 

V 

CBA-Arkansas 

2,452,988.82615 . 

100.0 

V 

CBA-Arkansas 

2,452,990.71092 . 

107.5 

V 

CBA-Arkansas 

2,452,990.83218 . 

108.0 

V 

CBA-Arkansas 

2,452,991.70918 . 

111.5 

V 

CBA-Arkansas 

2,452,991.83667 . 

112.0 

V 

CBA-Arkansas 

2,452,992.71317 . 

115.5 

V 

CBA-Arkansas 

2,452,992.83898 . 

116.0 

V 

CBA-Arkansas 

2,453,032.71822 . 

275.0 

B 

CBA-Arkansas 

2,453,033.72253 . 

279.0 

B 

CBA-Arkansas 

2,453,035.72843 . 

287.0 

B 

CBA-Arkansas 

2,453,691.86440 . 

2903.0 

none 

CBA-NM3 

2,453,691.98820 . 

2903.5 

none 

CBA-NM3 

2,453,693.87070 . 

2911.0 

none 

CBA-NM3 

2,453,693.99380 . 

2911.5 

none 

CBA-NM3 

2,453,694.87390 . 

2915.0 

none 

CBA-NM3 

2,453,694.99770 . 

2915.5 

none 

CBA-NM3 

2,453,695.75052 . 

2918.5 

none 

CBA-Utah 

2,453,695.87817 . 

2919.0 

none 

CBA-Utah 

2,453,696.00030 . 

2919.5 

none 

CBA-NM3 

2,453,696.00045 . 

2919.5 

none 

CBA-Utah 

2,453,696.88030 . 

2923.0 

none 

CBA-NM3 

2,453,697.75706 . 

2926.5 

none 

CBA-Utah 

2,453,697.88463 . 

2927.0 

none 

CBA-Utah 

2,453,698.00781 . 

2927.5 

none 

CBA-Utah 

2,453,700.89360 . 

2939.0 

none 

CBA-NM3 

2,453,701.01730 . 

2939.5 

none 

CBA-NM3 

2,453,702.89960 . 

2947.0 

none 

CBA-NM3 

2,453,703.77660 . 

2950.5 

none 

CBA-NM3 

2,453,703.90330 . 

2951.0 

none 

CBA-NM3 

2,453,704.77960 . 

2954.5 

none 

CBA-NM3 

2,453,704.90620 . 

2955.0 

none 

CBA-NM3 

2,453,705.77857 . 

2958.5 

none 

CBA-Utah 

2,453,705.78390 . 

2958.5 

none 

CBA-NM3 

2,453,705.90960 . 

2959.0 

none 

CBA-NM3 

2,453,705.91027 . 

2959.0 

none 

CBA-Utah 

2,453,706.78640 . 

2962.5 

none 

CBA-NM3 

2,453,706.91310 . 

2963.0 

none 

CBA-NM3 

2,453,708.79220 . 

2970.5 

none 

CBA-Utah 

2,453,708.92054 . 

2971.0 

none 

CBA-Utah 

2,453,709.79660 . 

2974.5 

none 

CBA-NM3 

2,453,709.92290 . 

2975.0 

none 

CBA-NM3 

2,453,710.79990 . 

2978.5 

none 

CBA-NM3 

2,453,710.92600 . 

2979.0 

none 

CBA-NM3 

2,453,711.80206 . 

2982.5 

none 

CBA-Utah 

2,453,711.80320 . 

2982.5 

none 

CBA-NM3 

2,453,711.92940 . 

2983.0 

none 

CBA-NM3 

2,453,711.93056 . 

2983.0 

none 

CBA-Utah 

2,453,712.80650 . 

2986.5 

none 

CBA-NM3 

2,453,712.93250 . 

2987.0 

none 

CBA-NM3 

2,453,713.80960 . 

2990.5 

none 

CBA-NM3 

2,453,713.93600 . 

2991.0 

none 

CBA-NM3 

TABLE  1 — Continued 


HJD 

Phase  Cycle 

Filter 

Observatory 

2,453,716.44484 . 

3001.0 

none 

CBA-Lesve 

2,453,719.83070 . 

3014.5 

none 

CBA-NM3 

2,453,719.95590 . 

3015.0 

none 

CBA-NM3 

2,453,725.59880 . 

3037.5 

none 

CBA-NM3 

2,453,725.72410 . 

3038.0 

none 

CBA-NM3 

2,453,725.84950 . 

3038.5 

none 

CBA-NM3 

2,453,725.97470 . 

3039.0 

none 

CBA-NM3 

2,453,729.61210 . 

3053.5 

none 

CBA-NM3 

2,453,729.73760 . 

3054.0 

none 

CBA-NM3 

2,453,729.86210 . 

3054.5 

none 

CBA-NM3 

2,453,729.98800 . 

3055.0 

none 

CBA-NM3 

2,453,730.61520 . 

3057.5 

none 

CBA-NM3 

2,453,730.74040 . 

3058.0 

none 

CBA-NM3 

2,453,730.86540 . 

3058.5 

none 

CBA-NM3 

2,453,730.99070 . 

3059.0 

none 

CBA-NM3 

2,453,731.61840 . 

3061.5 

none 

CBA-NM3 

2,453,733.50049 . 

3069.0 

none 

CBA-Lesve 

2,453,733.61898 . 

3069.5 

none 

CBA-Lesve 

2,453,740.39391 . 

3096.5 

none 

CBA-Lesve 

2,453,741.52569 . 

3101.0 

none 

CBA-Lesve 

2,453,742.52894 . 

3105.0 

none 

CBA-Lesve 

2,453,759.33366 . 

3172.0 

none 

CBA-Lesve 

2,453,677.82063 . 

2847.0 

V 

CBA-Arkansas 

2,453,677.95151 . 

2847.5 

V 

CBA-Arkansas 

2,453,679.951055 . 

2855.5 

V 

CBA-Arkansas 

2,453,680.82975 . 

2859.0 

B 

CBA-Arkansas 

2,453,680.95402 . 

2859.5 

V 

CBA-Arkansas 

2,453,681.83129 . 

2863.0 

B 

CBA-Arkansas 

2,453,681.83137 . 

2863.0 

V 

CBA-Arkansas 

2,453,684.84244 . 

2875.0 

B 

CBA-Arkansas 

2,453,684.84221 . 

2875.0 

V 

CBA-Arkansas 

2,453,684.96246 . 

2875.5 

V 

CBA-Arkansas 

2,453,685.84576 . 

2879.0 

B 

CBA-Arkansas 

2,453,685.96874 . 

2879.5 

V 

CBA-Arkansas 

2,453,685.84333 . 

2879.0 

V 

CBA-Arkansas 

2,453,686.84642 . 

2883.0 

V 

CBA-Arkansas 

Previously  published  in  Krajci  (2006). 


be  found  at  the  CBA  Web  site.4  These  light  curves  were  used  to 
obtain  measurements  of  the  times  of  minimum  for  both  primary 
and  secondary  eclipses  that  can  be  found  in  Table  1 . 

2.1.2.  USNO  Photometry 

Photometric  observations  obtained  at  the  US  Naval  Obser¬ 
vatory  (USNO),  Flagstalf  Station,  were  acquired  using  the  1  m 
Ritchey-Chretien  telescope.  The  telescope  was  equipped  with  a 
Tektronix  2048  x  2048  thinned  back-illuminated  CCD  and  Johnson 
UBV filters.  The  camera  has  a  scale  of  0.68"  pixel-1,  which  results 
in  a  field  of  view  of  23.2'  x23.2'. 

The  field  centered  on  FS  Aur  was  observed  on  three  consec¬ 
utive  nights,  2005  December  4-6  (UT),  in  the  B,  V,  and  U filters, 
respectively.  The  resulting  light  curves  are  shown  in  Figures  1  -3. 
Each  set  of  nightly  observations  was  preceded  with  a  set  of  20- 
50  bias  frames  and  1 0  dome  flats  in  all  three  filters.  A  total  of  1 1 0, 
176,  and  73  individual  exposures  were  obtained  in  the  B,  V,  and  U 
filters,  respectively.  Each  night  the  field  was  observed  during  an 
8  hr  window,  which  corresponded  to  hour  angles  between  — 4h  to 
+4h.  Although  each  night  was  devoted  to  a  single  filter,  at  least 
one  set  of  exposures  in  all  three  filters  was  also  acquired  each 
night.  The  typical  seeing  during  the  USNO  observing  run  was 
between  2"  and  3". 


4  See  http://cba.phys.columbia.edu. 
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HJD-2453700 

Fig.  1. —  {/-band  photometry  from  2005  December  6  obtained  with  the 
USNO  1  m  telescope.  The  differential  U  magnitudes  have  been  converted  to 
U  magnitudes  using  the  photometric  standard  stars  for  the  FS  Aur  field.  The 
residual  magnitudes  of  the  C  and  K  comparison  stars  are  also  shown. 


HJD-2453700 


Fig.  3. —  F-band  photometry  from  2005  December  5  obtained  with  the 
USNO  1  m  telescope.  The  differential  F  magnitudes  have  been  converted  to 
F  magnitudes  using  the  photometric  standard  stars  for  the  FS  Aur  field.  The 
residual  magnitudes  of  the  C  and  K  comparison  stars  are  also  shown. 


2.1.3.  R-Band  Photometry 

Time-series  differential  R-band  photometry  of  FS  Aur-79  was 
obtained  on  2006  January  2  at  the  Whispering  Pines  Observatory 
in  Harrison,  Arkansas,  using  a  0.4  m  f4.2  Newtonian  telescope 
with  a  Johnson  R  filter  and  SBIG  ST-6  CCD  camera.  The  differ¬ 
ential  R-band  light  curve  is  shown  in  Figure  4.  The  R  data  are  very 
valuable  as  the  extra  color  information  helps  reduce  the  unique¬ 
ness  problem  and  limit  parameter  space  in  modeling  binary  light 
curves,  especially  those  with  starspots. 

2.1.4.  Photometric  Reductions 

All  of  the  CCD  images  containing  photometry  data  were  pro¬ 
cessed  using  Cmuniwin.5  This  PC-based  software  with  a  graph¬ 
ical  user  interface  was  converted  from  Linux  C  programs  originally 
developed  by  Hroch  ( 1998).  In  short,  the  images  are  (1)  converted 
to  FITS  if  necessary,  (2)  flat-fielded  and  dark-subtracted  if  de¬ 
sired,  and  (3)  processed  to  find  stellar  targets  and  photometrically 
measured  using  the  algorithms  of  DAOPHOT  (Stetson  1987). 
Then,  (4)  target  lists  are  pattern  matched  to  identify  stars  in 


5  Available  at  http://integTal.sci.muni.cz/cmunipack/index.html. 


each  image  via  the  algorithm  of  Groth  (  1 986),  and  (5)  variable, 
comparison,  and  check  stars  are  selected  to  generate  differential 
photometry  and  light  curves. 

FS  Aur-79  lies  close  to  the  cataclysmic  variable  FS  Aur 
(Robertson  et  al.  2004),  and  the  field  has  been  observed  to  estab¬ 
lish  secondary  standard  stars  (Henden  &  Honeycutt  1997;  A.  A. 
Henden  2006,  private  communication).  Therefore,  there  are  com¬ 
parison  and  check  stars  in  the  field  with  standard  U,  B,  and  V 
magnitudes  but  not  R  as  of  yet.  The  magnitude  of  the  compari¬ 
son  star  was  added  to  the  differential  photometry  values  to  yield 
the  magnitudes  displayed  in  Figures  1-3. 

2.2.  Spectroscopy 
2.2.1.  WIYN 

Three  spectra  were  obtained  on  2006  February  22  UT  with  the 
HYDRA  Multiple  Object  Spectrograph  on  the  WIYN6  3.5  m 
telescope  at  Kitt  Peak  National  Observatory.  Table  2  lists  the 


6  The  WIYN  observatory  is  a  joint  facility  of  the  University  of  Wisconsin- 
Madison,  Indiana  University,  Yale  University,  and  the  National  Optical  Astron¬ 
omy  Observatory. 


Fig.  2. — 5-band  photometry  from  2005  December  4  obtained  with  the 
USNO  1  m  telescope.  The  differential  B  magnitudes  have  been  converted  to 
B  magnitudes  using  the  photometric  standard  stars  for  the  FS  Aur  field.  The 
residual  magnitudes  of  the  C  and  K  comparison  stars  are  also  shown. 


Fig.  4. — 5-band  photometry  from  2006  January  2  obtained  at  Whispering 
Pines  Observatory.  The  differential  R  magnitudes  are  shown,  since  comparison 
stars  have  not  been  calibrated  for  this  field.  The  residual  magnitudes  of  the  C 
and  K  comparison  stars  are  also  shown. 
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TABLE  2 

WIYN  Observations 


HJD 

Orbital  Phase 

2,453,788.6109458 . 

0.732 

2,453,788.6194817 . 

0.766 

2,453,788.6283672 . 

0.801 

HJD  and  orbital  phase  for  each  of  the  spectra  obtained.  The 
600  mm-1  grating  blazed  at  13.9°  was  used  in  the  first  order  to 
obtain  a  spectrum  covering  «535. 0-825.0  nm  pixel-1  with  a 
resolution  of  «0.14  nm  pixel-1.  Exposure  times  for  the  object 
spectra  were  600  s,  with  wispy  cirrus  clouds  present.  Dome-flat 
and  HeArNe  arc  lamp  images  were  obtained  for  calibrations  to 
remove  individual  CCD  pixel  variations  and  to  convert  pixels  to 
wavelengths.  Reductions  used  standard  IRAF7  onedspedtwodspec 
procedures  to  (1)  create  average  flat-field  and  bias  frames  ( im - 
combine),  (2)  subtract  the  bias  from  object  and  calibration  im¬ 
ages,  as  well  as  flat-field  correct  the  images  ( imarith ),  (3)  identify, 
trace,  sum,  and  extract  one-dimensional  apertures  from  the  two- 
dimensional  CCD  images  ( apall ),  (4)  determine  and  apply  the 
wavelength  calibration  to  object  apertures  from  the  HeArNe  lamp 
apertures  ( identify ,  refspec),  (5)  apply  corrections  to  eliminate  the 
motion  of  the  observatory  ( rvcorr ,  dopcor),  and  (6)  correct  the 
spectra  to  a  flat  continuum  for  relative  comparison  ( continuum ). 
Although  some  unassigned  fibers  yielded  sky,  the  sky  was  not  ex¬ 
tracted  or  subtracted. 

Spectra  of  FS  Aur-62  and  FS  Aur-63  were  also  obtained  si¬ 
multaneously  through  fibers,  as  these  are  the  comparison  stars 
used  in  the  differential  photometry.  The  photometric  calibrations 
of  these  stars  (Henden  &  Honeycutt  1997)  show  that  the  B  —  V 
colors  ofthese  stars  are  0.412  and  0.450,  respectively.  To  relative 
flux  calibrate  the  FS  Aur-79  spectra  we  used  the  spectra  of 
comparison  stars  with  similar  B  —  V  colors:  LTT  377  (B  —  V  = 
0.48),  LTT  1788  (B  -  V  =  0.47),  and  LTT  2415  (B  -  V  =  0.40). 
We  then  used  spectra  of  FS  Aur-62  and  FS  Aur-63  with  the  cal¬ 
ibrations  of  those  spectroscopic  standards  to  obtain  a  relative 
flux  calibration  of  the  FS  Aur-79  spectrum  using  IRAF  tasks 
standard,  sensfunc,  and  calibrate. 

The  average  spectrum  of  FS  Aur-79  is  shown  in  Figure  5  and 
corresponds  to  orbital  phase  «0.766.  The  lines  were  identified 
using  the  line  list  for  K  and  M  stars  from  Kirkpatrick  et  al.  ( 1991). 

2.2.2.  HET  HRS  Spectroscopy 

Radial  velocity  measurements  were  needed  to  determine  the 
individual  stellar  masses.  Given  the  approximately  6  hr  orbital 
period,  a  visual  magnitude  of  about  14,  and  the  need  to  acquire 
high-resolution  spectra,  a  large-aperture  telescope  was  needed. 
We  therefore  obtained  10  hr  of  priority-2  time  during  2005  se¬ 
mester  B  (September  through  January)  on  the  High  Resolution 
Spectrograph  ( HRS)  at  the  9.2  m  Hobby-Eberly  Telescope  (HET) 
through  NOAO.  Due  to  the  nature  of  the  HET  service  observing 
scheduling  only  one  object  spectrum  per  half  night  was  possible. 
In  order  to  get  as  complete  a  sampling  over  an  orbit  as  possible, 
we  requested  that  observations  be  done  as  randomly  in  time  as 
possible  over  the  semester,  avoiding  consecutive  nights  at  the 
same  time  of  night  due  to  the  approximately  6  hr  period  of  the  orbit. 
The  HRS  configuration,  yielding  a  resolving  power  of  15,000  and 
with  a  wavelength  coverage  of  approximately  4076-7838  A,  a 


7  IRAF  is  distributed  by  the  National  Optical  Astronomy  Observatory,  which 
is  operated  by  the  Association  of  Universities  for  Research  in  Astronomy,  Inc., 
under  cooperative  agreement  with  the  National  Science  Foundation. 


Fig.  5. — The  1.4  A  pixel-1  spectrum  of  FS  Aur-79  from  2006  February  21 
obtained  with  the  WIYN  telescope.  This  spectrum  was  obtained  at  an  orbital 
phase  of  about  0.76.  The  major  spectral  features  are  labeled  and  are  character¬ 
istic  of  K  and  M  stars. 

fiber  with  an  angular  diameter  on  the  sky  of  3",  one  sky  fiber,  no 
image  sheer,  no  iodine  gas  cell,  and  binning  of  2  x  5  for  3.2  pixels 
per  resolving  elements,  was  used.  Bias,  flat,  and  Th-Ar  compar¬ 
ison  lamp  frames  were  also  obtained  by  the  night  observers  for 
reduction  and  calibration  purposes.  Spectral  standards  were  also 
obtained,  giving  us  the  option  to  relative  flux  calibrate  the  spectra. 
Table  3  lists  the  details  of  the  24  HET  HRS  exposures. 

Standard  IRAF  tasks  were  used  to  process  the  frames  into  one¬ 
dimensional  spectra  normalized  to  the  continuum.  The  following 
reduction  steps  were  used  for  both  the  red  and  blue  channels: 
zerocombine,  flatcombine,  ccdproc  to  trim  and  apply  the  zero 


TABLE  3 

HET  Observations 


UT  Date 

UT  Time 

Orbital  Phase 

Exp 

(s) 

Skya 

S/Nb 

2005  Oct  9 . 

09:09 

0.609 

980 

S 

75 

2005  Oct  12 . 

08:31 

0.465 

900 

P 

124 

2005  Oct  20 . 

07:48 

0.246 

900 

P 

40 

2005  Nov  1 . 

07:29 

0.040 

900 

P 

78 

2005  Nov  3 . 

07:15 

0.976 

900 

P 

82 

2005  Nov  11 . 

06:40 

0.777 

900 

P 

150 

2005  Nov  11 . 

11:40 

0.609 

900 

P 

135 

2005  Nov  19 . 

06:13 

0.599 

900 

P 

140 

2005  Nov  19 . 

11:22 

0.454 

900 

P 

137 

2005  Nov  26 . 

05:40 

0.419 

900 

S 

131 

2005  Nov  26 . 

10:45 

0.263 

900 

P 

153 

2005  Dec  10 . 

10:01 

0.962 

900 

P 

120 

2005  Dec  19 . 

03:47 

0.810 

1000 

P 

135 

2005  Dec  19 . 

09:15 

0.717 

900 

P 

161 

2005  Dec  27 . 

03:36 

0.674 

900 

s 

114 

2005  Dec  28 . 

09:06 

0.574 

900 

s 

93 

2006  Jan  3 . 

02:58 

0.478 

900 

p 

114 

2006  Jan  11 . 

02:21 

0.269 

900 

s 

110 

2006  Jan  11 . 

07:42 

0.158 

900 

s 

130 

2006  Jan  19 . 

07:34 

0.030 

900 

s 

2006  Jan  20 . 

01:51 

0.067 

900 

p 

90 

2006  Jan  21 . 

07:09 

0.935 

900 

p 

2006  Jan  21 . 

07:27 

0.984 

900 

p 

72 

2006  Jan  29 . 

01:35 

0.903 

900 

p 

89 

a  S:  spectroscopic;  P:  photometric. 
b  Per  resolution  element. 
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Fig.  6. —  Spectra  of  the  Ba  i — Ti  i— Fe  i-Ca  i  blend  at  representative  orbital  phases  obtained  with  the  HET  FIRS  (order  94)  normalized  to  the  continuum  ( horizontal 
dotted  lines).  The  velocity  scale  is  for  the  center-of-mass  frame  of  the  binary. 


correction,  cosmicrays ,  ctpall  to  find  and  define  the  apertures  for 
each  spectral  order  on  the  red  and  blue  chips,  apscatter  to  remove 
any  scattered  light,  apflatten  to  get  a  flat-field  solution  using  the 
combined  flat,  ccdproc  to  apply  the  flattening  to  the  spectral  frames, 
and  apall  once  again  to  extract  the  spectral  orders.  HET  HRS  orders 
78-145  were  extracted,  corresponding  to  spectral  coverage  from 
4200  to  7900  A.  When  necessary,  sky  subtraction  was  done  for 
those  nights  near  the  full  moon  when  the  sky  background  was 
significant  and  for  those  orders  that  have  problematic  sky  lines 


(i.e.,  order  104,  Na  i  D  doublet).  Wavelength  calibration  was 
done  for  each  spectral  order  of  interest  using  ident,  reident, 
refspec,  and  dispcor.  The  spectra  were  then  converted  to  the  helio¬ 
centric  frame  of  reference  using  rvcor  and  dopcor.  The  instru¬ 
mental  flux  was  then  normalized  to  the  continuum  by  using  either 
contin,  if  it  could  correctly  find  the  continuum,  or  by  doing  the 
continuum  fitting  interactively  using  the  option  under  splot.  Spec¬ 
tra  from  order  94  are  shown  in  Figure  6,  sorted  by  phase.  The 
main  spectral  feature  in  this  order  is  the  absorption  line  blend  of 


Velocity  (km/sec) 


-300-200-100  0  100  200  300 


Wavelength  [A] 

Fig.  7. —  Spectra  oftheNai  D  doublet  at  representative  orbital  phases  obtained  with  the  HET  HRS  (order  104)  normalized  to  the  continuum  (horizontal  dotted  lines). 
The  velocity  scale  is  for  the  center-of-mass  frame  of  the  binary. 


No.  5,  2007 


ECLIPSING  BINARY  [HH97]  FS  AUR-79 


1939 


0.006 


0.004 


0.002 


O 

I  0 
O 


-0.002 


-0.004 


-0.006 


Cycle 


Fig.  8. — Difference  between  observed  and  calculated  (O  —  C)  eclipse  tim¬ 
ings  from  CCD  photometry  of  FS  Aur-79  obtained  by  members  of  the  CBA. 
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Fig.  10. — Phased  5-band  photometry  from  2005  December  4  obtained  with 
the  USNO  1  m  telescope.  Solid  curve :  Best-fit  Nightfall  model.  Dashed  curve : 
BM3  model. 

eris,  min  I  =  HJD  2,452,963.74445  +  0.250816A,  are  shown  in 
Figure  8.  The  fact  that  (O  —  C)  has  been  approximately  zero  over 
four  years  of  photometry  would  indicate  no  changes  in  the  orbital 
period  and  therefore  no  indication  of  any  significant  mass  transfer. 


Ti  1 26497.689,  Ba  1 26498.759,  Fe  1 26498.950,  Ca  1 26499.694, 
and  Fe  1 26501.681.  Figure  7  shows  spectra  from  order  104  con¬ 
taining  the  Na  i  D  doublet  sorted  by  phase. 

3.  ANALYSIS 

3.1.  Ephemeris 

Both  primary  and  secondary  eclipse  times  of  minimum  were 
determined  using  the  method  ofKwee  &  van  Woerden  (1956) 
and  are  listed  in  Table  1 .  All  of  these  timings  were  then  used  in  a 
linear  least-squares  fit  of  the  form  T  =  Tq  +  PN  to  refine  the  pe¬ 
riod  and  ephemeris,  where  Tq  is  the  Heliocentric  Julian  Date 
of  primary  eclipse  minimum  (occultation  of  the  more  massive, 
hotter  star)  from  Robertson  et  al.  (2004),  N  is  the  orbital  cycle 
since  the  epoch  T0,  P  is  the  orbital  period,  and  T  values  are  the 
predicted  times  of  primary  eclipses.  This  yielded  a  period  of 
0.250816(1)  days.  The  differences  between  observed  eclipses 
and  their  calculated  timings  (O  —  C )  based  on  this  new  ephem- 


3.2.  Light  Curves 

The  U,  B,  V,  and  R  phased  light  curves  of  FS  Aur-79  are  dis¬ 
played  in  Figures  9-12.  The  morphology  of  these  light  curves 
place  FS  Aur-79  in  the  category  of  (3  Lyrae  eclipsing  binaries,  since 
the  light  curves  vary  continuously  and  the  primary  (phase  =  0.00) 
and  secondary  eclipses  (phase  =  0.50)  are  of  unequal  depths. 
Table  4  lists  the  U,  B,  V,  and  R  magnitudes  at  primary  eclipse, 
secondary  eclipse,  max  I  (phase  =  0.25),  and  max  II  (phase  = 
0.75).  By  comparing  the  magnitude  of  the  secondary  eclipse  to 
the  magnitudes  during  the  quadratures  we  are  able  to  find  a  lower 
limit  to  the  magnitude  differences  between  the  primary  and  sec¬ 
ondary  star,  since  we  do  not  know  if  the  eclipses  are  total.  The 
primary  star  is  at  least  1 .42  mag  brighter  than  the  secondary  in  V 
and  at  least  1.33  mag  brighter  than  the  secondary  in  B. 

The  light  curves  exhibit  the  O’Connell  effect  (unequal  max¬ 
ima)  of  BY  Dra  variables,  indicating  starspots.  The  amount  of  the 
O’Connell  effect  in  each  color  can  be  quantified  by  the  difference 


Fig.  9. — Phased  £/-band  photometry  from  2005  December  6  obtained  with  Fig.  11. —  Phased  F-band  photometry  from  2005  December  5  obtained  with 

the  USNO  1  m  telescope.  Solid  curve :  Best-fit  Nightfall  model.  Dashed  curve :  the  USNO  1  m  telescope.  Solid  curve :  Best-fit  Nightfall  model.  Dashed  curve : 

BM3  model.  BM3  model. 


1940 


AUSTIN  ET  AL. 


Vol.  133 


Lightcurve  Phase 

Fig.  12. — Phased  i?-band  photometry  from  2006  January  2  obtained  at 
Whispering  Pines  Observatory.  Solid  curve :  Best-fit  Nightfall  model.  Dashed 
curve :  is  the  BM3  model. 


(Am)  between  max  II  and  max  I  and  are  listed  in  Table  5.  Also 
shown  in  Table  5  are  the  Am  U  —  B,  B  —  V ,  and  V  —  R  colors. 

Color  curves  were  computed  by  first  interpolating  each  light 
curve  to  a  phase  grid  with  phase  increments  of  0.0125,  then  sub¬ 
tracting  the  light  curves  of  different  colors,  and  therefore  also 
different  nights.  Figure  13  shows  the  U  —  B,  B  —  V,  and  dif¬ 
ferential  V  —  (R  —  comp)  colors  with  respect  to  orbital  phase. 
Table  4  includes  the  U  —  B  and  B  —  V  colors  at  primary  eclipse, 
secondary  eclipse,  and  max  I  and  max  II.  During  secondary 
eclipse  the  secondary  star  is  mostly  hidden,  and  the  observed 
color  will  be  for  the  primary  star.  B  —  V  during  secondary  eclipse 
is  about  1 .29,  which  would  correspond  to  a  main-sequence  star 
with  an  effective  temperature  of  around  4100  K  and  spectral  type 
of  about  K7.  During  secondary  eclipse,  JJ  —  B  =  0.860,  which 
would  correspond  to  a  main-sequence  star  with  effective  tem¬ 
perature  of  about  4700  K.  This  U  —  B  excess  is  an  indicator  of 
chromospheric  activity  (which  is  also  seen  in  our  spectra),  and 
therefore  this  U  —  B  cannot  be  used  to  determine  an  accurate 
photospheric  temperature  (Amado  2003). 

The  apparent  V  magnitude  of  the  primary  is  13.829,  and  the 
B  —  V  color  is  1.292.  Assuming  the  primary  is  a  main-sequence 
star,  then  the  distance  modulus  V  —  M  =  5.82,  and  the  distance 
is  about  146  pc  (similar  to  the  distance  to  the  Pleiades).  Given 
this  distance,  reddening  and  extinction  are  probably  small  enough 
to  be  ignored. 

3.3.  Spectra 

3.3.1.  Spectral  Classification 

The  WIYN  spectrum  was  obtained  when  the  orbital  phase 
of  the  system  was  about  0.76.  Figure  14  shows  our  WIYN  spec¬ 
tra,  a  main-sequence  K7  spectrum  (star  FID  237903),  and  a  main- 
sequence  Ml  spectrum  (star  HD  204445)  from  the  Indo-US 


Library  of  Coude  Feed  Spectra8  (Valdes  et  al.  2004).  The  K7  and 
Ml  spectra  were  broadened  with  a  rotational  kernel  calculated 
for  a  v  sin  i  value  of  1 35  km  s  1  (approximate  rotational  speed  of 
FS  Aur-79)  and  a  limb-darkening  coefficient  of  0.6  (see,  e.g., 
Gray  2005,  p.  465).  The  match  between  the  WIYN  spectra  and 
the  K7  V  spectrum  is  very  close  between  5200  and  7000  A.  Be¬ 
tween  7000  and  7500  A  the  WIYN  spectrum  looks  a  bit  more 
like  an  early  dM  spectrum.  But  the  spectrum  between  7000  and 
7500  A  is  mostly  a  blend  of  a  K7  V  spectrum  and  an  M-type 
spectrum.  Given  the  uncertainties  involved  in  the  atmospheres  of 
cool  M-type  stars  and  the  lack  of  absolute-flux  calibration  avail¬ 
able,  it  is  difficult  at  this  time  to  accurately  pinpoint  the  precise 
spectral  type  of  the  secondary  from  the  WIYN  spectrum.  Binary 
modeling  in  §  3.3  indicates  a  secondary  temperature  consistent 
with  an  M4  V  star,  However,  the  mass  and  radius  found  by  our 
binary  modeling  are  more  consistent  with  approximately  an 
M2  V  or  M3  V  star  by  comparing  to  the  masses  and  radii  de¬ 
termined  for  YY  Gem  and  CM  Dra  (Strassmeier  et  al.  1 988).  Of 
all  the  potential  spectroscopic  chromospheric  indicators  that 
are  accessible  in  our  wavelength  coverage,  such  as  emission  in 
hydrogen  Balmer  lines,  Ca  n  H,  K  emission  cores,  He  i  25876 
absorption  or  emission  core,  and  Na  i  D  doublet  emission  cores, 
this  star  shows  strong  emission  in  Ha  along  with  correspond¬ 
ing  emission  in  H /3.  We  therefore  estimate  the  spectral  types  as 
dK7e  for  the  primary  star  and  dM3e  for  the  secondary  star. 

3.3.2.  Velocities 

FS  Aur-79  is  a  single-lined  spectroscopic  binary,  with  the  lines 
from  the  primary  star  most  obvious,  and  the  presence  of  the  lines 
of  the  secondary  detected  by  their  influence  on  the  line  profiles  and 
TiO  band  strengths.  The  luminosity  difference  (in  the  V  band  the 
secondary  is  at  most  0.27  times  as  bright  as  the  secondary),  the 
short  period,  and  the  short  exposure  times  prohibited  us  from  re¬ 
solving  the  lines  of  the  secondary. 

The  7-velocity  (radial  velocity  of  the  center  of  mass  of  the 
system  relative  to  the  Sun)  and  the  radial  velocities  of  the  primary 
can  be  measured  directly.  We  first  determined  the  7-velocity 
by  choosing  lines  that  we  were  confident  were  not  blends,  for 
example,  Ca  1  lines  whose  rest  wavelengths  are  6122.219  and 
6162.172  A.  Wavelengths  of  spectral  features  of  interest  were 
measured  using  the  equivalent  width  function  in  the  1RAF  task 
splot.  A  sine  wave  was  fit  to  the  wavelength  versus  orbital  phase 
for  these  lines  that  determined  the  wavelength  of  the  line  in  the 
center-of-mass  frame  relative  to  the  Sun,  6123.94  ±  0.07  and 
6163.78  ±0.11  A,  respectively,  resulting  in  7-velocities  of 
+84.3  ±  3.6and+78.3  ±  5.4  kms-1.  The  chromospheric  Ha  line 
does  not  show  any  modulation  of  its  wavelength  as  a  function  of 
orbital  phase,  so  it  too  was  used  to  determine  the  7-velocity.  Here 
a  line  with  slope  zero  was  fit  to  the  wavelength  versus  orbital 
phase  data,  resulting  in  an  observed  wavelength  of  6564.68  ± 
0.12  A  and  a  7-velocity  of +85. 2  ±5.5  kms-1.  We  conclude  that 


8  See  http://www.noao.edu/cflib/. 


TABLE  4 


Magnitudes  and  Colors  during  Eclipses  and  Quadratures 


Phase 

u 

B 

V 

R  —  comp 

U-B 

B  —  V 

0.00  (primary  eclipse) . 

16.579  +  0.007 

15.811+0.007 

14.479  +  0.06 

0.2183  +  0.015 

0.768 

1.33 

0.25  (max  I) . 

15.739  +  0.007 

14.892  +  0.007 

13.620  +  0.06 

-0.61155  +  0.015 

0.847 

1.272 

0.50  (secondary  eclipse) . 

15.970  +  0.007 

15.121+0.007 

13.829  +  0.06 

-0.34813  +  0.015 

0.849 

1.292 

0.75  (max  II) . 

15.595  +  0.007 

14.797  +  0.007 

13.517  +  0.06 

-0.68721  +0.015 

0.798 

1.280 
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TABLE  5 

O’Connell  Effect 


Quantity  Max  II  —  Max  I  Color 


Amu .  —0.144 

A% .  —0.095 

Amy .  —0.103 

A  mR .  -0.0757 

Amu  —  Amg .  •••  —0.049 

A  mg  —  Amy .  . . .  0.008 

Amy  —  Amu .  . . .  —0.027 


the  7-velocity  of  the  system,  based  on  the  mean  of  the  values 
above,  is  +82.6  ±5.5  km  s_1. 

We  ran  the  1RAF  cross-correlation  task  fxcor  on  HET  HRS 
orders  94  through  100,  covering  a  wavelength  range  of  6070- 
6535  A,  which  contains  several  strong  lines.  We  used  a  good  S/N 
spectrum  near  orbital  phase  0.5  as  the  template.  The  resulting 
radial  velocity  curve  had  unacceptably  large  radial  velocity  er¬ 
rors.  The  large  errors  are  likely  the  result  of  the  spectra  changing 
with  phase  due  to  not  only  Doppler  shifting  but  also  due  to  star- 
spot  groups,  and  to  a  lesser  extent  the  changing  relative  con¬ 
tributions  to  the  light  from  the  two  stars  due  to  the  eclipses.  By 
measuring  the  wavelengths  of  over  a  dozen  strong  absorption 
lines  using  the  equivalent  width  function  in  the  1RAF  task  splot, 
we  were  not  as  sensitive  to  these  other  effects.  Figure  15  shows 
an  example  radial  velocity  curve  from  the  absorption  line  of  the 
blend  at  6500  A  in  the  center-of-mass  frame  of  the  system.  The 
error  bars  in  phase  indicate  the  exposure  time  for  each  spectrum, 
and  the  error  bars  in  velocity  indicate  the  spectral  resolution  of 
the  HRS  spectra.  The  maximum  radial  velocity  ( v  sin  i )  of  the  pri¬ 
mary  is  about  125  ±  15  km  s_1.  Also  shown  in  Figure  15  are  the 
radial  velocity  curves  of  both  the  primary  and  secondary  stars 
from  the  best-fit  binary  model  discussed  in  §  3.4  below. 

3.4.  Binary  Modeling 
3.4.1.  Optimization 

Light  and  radial  velocity  curves  were  fit  simultaneously  using 
the  program  Nightfall.9  Nightfall  has  built-in  optimization  pro- 


Available  at  http://www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall 
.html. 


Fig.  13. — Phased  U  —  B,  B  —  V,  and  V  —  (R  —  comp)  light  curve.  Solid 
lines  represent  the  best-fit  Nightfall  model. 


Fig.  14. — The  1.4  Apixel  1  spectrum  of  FS  Aur-79  from  2006  February  21 
obtained  with  the  WIYN  telescope  compared  to  rotationally  broadened  (v  sin  i  — 
135  km  s  _1)  spectra  of  K7  V and  Ml  V  stars  (Indo-US  Fibrary  stars  HD  237903 
and  HD  204445  [Valdes  et  al.  2004]). 

cedures  (simplex  and  simulated  annealing).  The  temperature  of 
the  primary  derived  from  the  photometry  and  spectroscopy  was 
set  to  4100  K.  The  orbital  period  was  set  to  0.2508  days.  Model 
atmosphere,  detailed  reflection,  and  a  square-root  limb-darkening 
law  options  were  enabled.  The  temperature  of  the  secondary  was 
initially  set  to  a  value  consistent  with  a  early  dM  star.  We  then  had 
Nightfall  search  parameter  space,  allowing  it  to  vary  the  temper¬ 
ature  of  the  secondary,  mass  ratio,  inclination  angle,  and  filling 
factors  for  the  primary  and  secondary.  Initially  only  spotless  mod¬ 
els  were  computed.  Nightfall  converged  on  the  values  of  mass 
ratio,  inclination  angle,  and  filling  factor  listed  in  Table  6.  Night¬ 
fall  computes  a  reduced  x2  (x2)  °f  the  simultaneous  fit  of  light 
and  radial  velocity  curves  for  the  best  spotless  model  of  approx¬ 
imately  54. 

Spotless  models  could  not  reproduce  the  O’Connell  effect,  and 
therefore  various  starspot  models  were  tried.  Starspots  are  to  be 
expected  for  stars  that  have  a  rotation  period  of  only  0.2508  days 
( Bopp  &  Fekel  1977;  Young  et  al.  1987).  The  temperature  of  the 
secondary  and  the  longitude,  size,  and  temperature  factor  relative 
to  the  photosphere  (dimming  factor)  of  the  circular  spot  were 
allowed  to  vary.  The  longitude  convention  used  in  Nightfall  has 
0°  on  the  secondary  facing  LI  and  0°  on  the  primary  facing  away 
from  LI .  An  optimized  model  with  a  single  equatorial  spot  on  the 


Fig.  15. — Radial  velocity  curve  of  the  primary  star.  Solid  curves  are  the 
model  radial  velocity  curves  for  the  primary  and  secondary. 
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TABLE  6 

Nightfall  Models 


Parameter 

Spotless 

One  Spot 

Three  Spot 

Uncertainty 

Stellar  Parameters 

Period  (days)3 . 

.  0.251 

0.251 

0.251 

i  (deg) . 

.  83.3 

83.3 

83.3 

±2.3 

q  (=M2/M,) . 

.  0.527 

0.527 

0.527 

±0.030 

Total  mass  (M0)a . 

.  0.9 

0.9 

0.9 

±0.1 

T\  (K)a . 

.  4100 

4100 

4100 

±25 

«  K » . 

.  3170 

3262 

3425 

±25 

Fill  factor  i . 

.  1.00 

1.00 

1.00 

±0.01 

Fill  factor 2 . 

.  0.950 

0.950 

0.950 

±0.030 

Gravity  darkening  i . 

.  0.082 

0.082 

0.082 

Gravity  darkening  2 . 

.  -0.016 

-0.002 

0.020 

Albedo  1  =  Albedo  2 . 

.  0.5 

0.5 

0.5 

Spot(s)  Primary  (Nightfall  Coordinates) 

Longitude- 1  (deg) . 

272 

310 

Latitude- 1  (deg) . 

0 

0 

Radius- 1  (deg) . 

66 

55 

±10 

Dimfactor-1 . 

0.985 

0.992 

±0.007 

Longitude-2  (deg) . 

210 

Latitude-2  (deg) . 

0 

Radius-2  (deg) . 

68 

±10 

Dimfactor-2 . 

0.971 

±0.010 

Spot(s)  Secondary  (Nightfall  Coordinates) 

Longitude . 

0 

Latitude . 

0 

Radius  (deg) . 

45 

±5 

Dimfactor . 

0.893 

±0.050 

Output  Parameters 

xl . 

.  54 

22 

9 

Mean  temperature  1  ( K) . 

.  4106 

4089 

4062 

±25 

Mean  temperature  2  (K) . 

.  3170 

3262 

3365 

±25 

Mean  radius  1  (Rq) . 

.  0.673 

0.673 

0.673 

Mean  radius  2  (Rq) . 

.  0.477 

0.477 

0.477 

Distance  (Rq) . 

.  1.616 

1.616 

1.616 

Orbital  velocity]  (km  s-1).... 

.  111.61 

111.61 

111.61 

Orbital  Velocity2  (km  s_1)  ... 

.  211.83 

211.83 

211.83 

Mi  (MX . 

.  0.59 

0.59 

0.59 

±0.02 

m2(MX . 

.  0.31 

0.31 

0.31 

±0.09 

a  Held  fixed. 


secondary  resulted  in  a  xl  °f  about  22.  The  best  model  param¬ 
eters  for  a  single  spot  on  the  secondary  are  also  listed  in  Table  6. 

Although  this  single-spot  model  is  an  improvement  over 
the  spotless  model,  one  spot  is  inadequate  to  account  for  the 
O’Connell  efFect  seen.  In  other  words,  the  photospheric  tempera¬ 
ture  inhomogeneity  is  more  complex  than  a  single  starspot.  In 
fact,  what  we  are  modeling  as  single  spots  are  probably  large  spot 
groups.  The  UB  VR  light  curves  were  fit  the  best  by  a  three-spot 
model  (Figs.  9-12),  one  circular  spot  on  the  secondary  and  two 
circular  spots  on  the  primary  (Fig.  16).  Again,  the  optimization 
varied  the  temperature  of  the  secondary  and  the  longitude,  size, 
and  dimming  factor  for  the  spots  and  resulted  in  a  fit  with  a  Xl  of 
about  9.  Table  6  also  shows  the  best-fit  parameters  for  the  three- 
spot  model  obtained  using  Nightfall.  The  uncertainties  listed  in 
Table  6  represent  the  parameter  ranges  over  which  the  minimum 
xl  remains  essentially  unchanged.  We  also  compared  the  mod¬ 
elled  and  observed  U  —  B,  B  —  V,  and  V  —  (R  —  comp)  color 
curves  (Fig.  13).  The  B  —  V  curve  fits  the  observed  values  closely 


at  all  phases.  The  U  —  B  curve  fits  generally  well,  but  the  observed 
data  show  features  that  are  not  reproduced  by  the  model.  The 
V  —  (R  —  comp)  also  fits  quite  well,  but  the  model  seems  to  have 
systematic  differences  compared  to  the  data.  The  actual  spots  are 
probably  spot  groups  that  are  irregular  in  shape,  and  therefore  the 
model  cannot  reproduce  in  detail  the  sunspot  group  shapes. 

The  model  parameters  optimized  by  Nightfall  were  cross¬ 
checked  by  modeling  with  BinaryMaker,  version  3.0  (BM3; 
Bradstreet  &  Steelman  2002).  BM3  assumes  the  stars  to  be  black- 
bodies,  whereas  the  stellar  atmosphere  option  was  used  for  the 
Nightfall  models.  To  achieve  a  fit  similar  to  that  obtained  with 
Nightfall  (Figs.  9-12)  we  had  to  change  the  inclination  angle 
slightly  and  use  a  secondary  temperature  50  K  cooler  when  us¬ 
ing  BM3;  otherwise,  the  same  parameter  values  were  used.  The 
equivalent  best-fit  three-spot  parameters  used  in  BinaryMaker 
are  shown  in  Table  7.  As  can  be  seen  in  Figures  9-12,  the  U and 
B  light  curves  from  BM3  do  not  fit  as  well  as  the  Nightfall  equiv¬ 
alent,  but  the  hand  R  light  curves  are  closer  to  being  identical.  We 
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Fig.  16. — Configuration  of  spots  for  the  best-fit  three-spot  binary  model  shown  at  orbital  phases  0.00,  0.25,  0.50,  and  0.75. 


TABLE  7 

BinaryMaker  Model 


Parameter  Three  Spot 


Stellar  Parameters 

Mass  ratio .  0.527 

Fillout  1 .  0.0 

Fillout  2 .  -0.01 

Temperature  1 .  4100 

Temperature  2 .  3375 

G1  =  G2 .  0.32 

XI .  0.887 

X2 .  0.973 

Reflection  1  =  Reflection  2 .  0.5 

Inclination .  81 

Period .  0.251 


Spots  Primary  (BinaryMaker  Coordinates) 

Longitude .  330 

Latitude .  90 

Spot  radius .  68.2 

Temperature  factor .  0.97 

Longitude .  230 

Latitude .  90 

Spot  radius .  54 

Temperature  factor .  0.99 


Spot  Secondary  (BinaryMaker  Coordinates) 

Longitude .  0.0 

Latitude .  90 

Spot  radius .  45 

Temperature  factor .  0.89 


attribute  these  differences  to  the  fact  that  model  atmospheres  were 
used  for  the  Nightfall  modeling  and  blackbodies  for  the  BM3 
model. 

By  simultaneously  fitting  the  U,  B,  V,  and  R  light  curves  we 
improved  the  likelihood  of  finding  a  unique  solution.  One  should 
interpret  the  starspots  in  the  model  simply  as  a  way  to  introduce 
temperature  differences  on  the  surfaces  of  the  stars,  since  one 
cannot  distinguish  a  cool  spot  from  a  hot  spot  on  the  opposite 
side  of  the  system.  By  comparing  the  observed  and  modelled 
color  curves  (U  —  B,  B  —  V,  and  V  —  R)  in  Figure  13,  one  can 
see  that  the  model  does  a  good  job  of  reproducing  the  observed 
temperature  inhomogeneities. 

3.4.2.  Spectral  Synthesis 

Given  the  spectral  types  of  the  component  stars  and  their 
observed  and  predicted  radial  velocities  we  have  synthesized  the 
Na  D  line  profile  for  zero  and  quadrature  orbital  phases.  This  pro¬ 
vides  a  qualitative  check  on  the  orbital  velocity  of  the  secondary 
star,  whose  spectral  features  are  not  noticeable  from  visual  in¬ 
spection  of  the  HET  HRS  spectra.  Figure  17  shows  the  observed 
and  synthesized  profiles  for  order  100.  A  K7  V  spectrum  and  an 
M 1  V  spectrum  were  normalized  to  the  continuum  over  the  wave¬ 
length  range  of  the  HET  HRS  order.  The  M 1  V  spectrum  was  then 
scaled  by  a  factor  of  about  0.225  (relative  luminosity  of  an  Ml  V 
star  in  units  of  a  K7  V  star’s  luminosity).  The  composite  spectrum 
was  then  normalized  to  the  continuum.  The  result  mimics  the  line 
profiles  of  the  sodium  doublet,  which  indicates  that  the  second¬ 
ary  star  spectrum  is  having  an  observable  effect  on  this  spectral 
feature  and  that  the  predicted  secondary  velocity  is  probably 
close  to  the  actual  velocity. 

4.  DISCUSSION 

FS  Aur-79  represents  a  very  near  contact  system  containing 
late-type  (K-M  dwarf)  components.  Near-contact  systems  have 
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Fig.  17. —  Sodium  doublet  and  synthesized  profiles  for  phases  0.0,  0.25,  0.5,  and  0.75. 


been  defined  by  Shaw  ( 1990),  in  which  the  binaries  have  a  pe¬ 
riod  of  less  than  1  day,  show  strong  tidal  interactions,  and  have 
surfaces  very  near  one  another  but  not  in  contact.  These  stars  are 
related  to  W  UMa  type  systems  and  are  probable  precursors  to 
such  systems.  W  UMa  systems  have  periods  less  than  1  day,  and 
their  light  curves  show  distortions  due  to  their  close  proximity, 
similar  to  near-contact  systems.  There  are  no  known  W  UMa 
type  systems  with  an  M  star  component,  so  this  system  repre¬ 
sents  a  unique  precursor  to  such  a  low-mass  W  UMa  system. 

Unlike  the  W  UMa  systems,  however,  the  near-contact  sys¬ 
tems  show  unequal  eclipse  depths  in  their  primary  and  secondary 
minima.  This  can  be  explained  with  the  fact  that  the  stars  have 
different  temperatures  and  are  therefore  not  in  contact  and  hence 
not  in  thermal  equilibrium,  indicative  of  / 3  Lyrae  systems. 

There  has  been  some  speculation  as  to  the  actual  existence  of 
W  UMa-type  systems  that  are  currently  not  in  the  contact  state 
due  to  thermal  relaxations  (Lucy  1976;  Flannery  1976;  Wang 
1994)  or  that  the  stars  are  somehow  able  to  maintain  temperature 
differences  within  a  contact  discontinuity  (Shu  et  al.  1976;  Wang 
1999).  These  proposed  systems  have  been  referred  to  as  the 
B-type  (apart  from  W-type  and  A-type)  subclass  of  W  UMa  sys¬ 
tems.  Indeed,  it  has  been  shown  that  some  of  the  observational 
candidates  for  the  B-type  W  UMa  systems  can  be  the  result  of 
spurious  effects  produced  when  starspots  are  not  included  in 
light-curve  models.  Maceroni  &  van’t  Veer  ( 1990)  showed  the 
effect  of  modeling  a  near-contact  system  with  and  without  starspots. 
Their  results  indicated  that  a  spotless  contact  model  could  repro¬ 
duce  equally  well  the  light  curves  of  spotted  noncontact  systems. 

Starspots  are  notorious  for  introducing  ambiguity  in  light- 
curve  modeling  in  chromospherically  active  binaries  such  as  RS 
C  Vn  and  Algols,  as  well  as  BY  Dra  rotational  variables,  especially 
when  there  is  only  one  color-band  light  curve.  Multiwavelength 
studies,  coupled  with  suitable  spectroscopy  (i.e.,  radial  velocities 
and  phase-correlated  chromospheric  activity  indicators)  can  often 
disentangle  multiple  light-curve  models  with  the  larger  number  of 
free  parameters  introduced  by  inclusion  of  starspots. 

Close  binaries  with  late-type  components  are  very  rare  for  two 
reasons.  The  first  and  most  obvious  is  their  intrinsic  faintness 


and,  therefore,  observational  bias  in  detecting  and  observing  them 
in  the  first  place.  The  other  is  due  to  their  evolutionary  history. 
Theories  of  protostellar  collapse  involving  formation  of  close  bi¬ 
nary  stars  with  separations  1  AU  or  less  are  problematic,  as  re¬ 
viewed  by  Bonnell(  2001).  Separations  of  a  few  AU  are  possible, 
however,  and  further  evolution  of  the  binary  can  lead  to  short 
orbital  period  systems  through  the  loss  of  angular  momentum. 
Therefore,  it  takes  time  for  these  systems  to  form. 

It  is  widely  known  that  young  stars  are  rapid  rotators  and  that 
main-sequence  evolution  causes  a  slowing  of  rotation.  The  de¬ 
crease  of  Ca  ii  H  and  K  emission  as  a  function  of  age  in  main- 
sequence  spectral  types  was  clearly  demonstrated  by  Wilson 
( 1963).  This  leads  to  an  activity  versus  age  relationship  among 
main-sequence  stars  as  their  stellar  winds,  which  are  coupled  to 
the  magnetic  field  and  are  forced  to  corotate  with  stellar  rotation 
out  to  large  radii,  extract  angular  momentum.  The  result  is  that 
the  stars  rotate  more  slowly  over  time.  This  in  turn  reduces  the 
stellar  activity  level,  as  the  speed  of  rotation  is  tied  to  magnetic 
dynamo  strength.  Since  the  result  of  extracting  angular  mo¬ 
mentum  from  a  binary  system  is  to  tighten  the  binary  orbit,  this 
angular  momentum  loss  mechanism,  when  operating  in  a  binary 
system,  results  in  the  evolution  of  the  binary  to  shorter  orbital 
periods.  Stepien  (1995)  has  formulated  the  angular  momen¬ 
tum  loss  via  a  magnetized  wind  and  the  spin-up  rate  of  binaries 
through  their  decreasing  orbital  periods  calibrated  from  observa¬ 
tions  of  the  spin-down  rates  of  single  stars  (Simon  et  al.  1985). 
These  results  show  that  the  timescales  for  angular  momentum 
loss  of  rapidly  rotating  single  stars  increase  substantially  for 
cooler  stars  and  later  spectral  types.  Therefore,  the  timescale 
for  contact  binary  formation  from  initially  more  widely  sepa¬ 
rated  components  will  increase  substantially  with  later  spectral  type 
components.  Even  the  closest  in  situ  formed  binaries  (2-5  AU) 
with  low  enough  masses  would  not  have  had  enough  time  to  reach 
the  near-contact  or  contact  stages  within  the  age  of  our  galaxy. 
Stepien’s  (1995)  results  show  that  for  a  binary  with  two  0.6  M0 
components,  angular  momentum  loss  operating  over  2.5  Gyr  is 
needed  to  form  a  contact  binary  having  an  initial  orbital  period  of 
2  days.  For  an  initial  orbital  period  of  3.4  days  this  increases  to 
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Fig.  18. — Examples  of  Hct  and  H (3  in  emission  lines  at  phases  0.0  and  0.5. 


Fig.  19. — TiO  bands  and  band  depths. 


6  Gyr.  The  primary  star  of  FS  Aur-79  has  a  mass  of  approxi¬ 
mately  0.6  M0  based  on  its  color  and  spectral  type.  Our  model 
indicates  a  mass  ratio  q  of  about  0.527,  which  results  in  a  sec¬ 
ondary  mass  of  approximately  0.3  M0,  which  is  consistent  with 
a  main-sequence  star  with  the  modelled  temperature  of  3425  K. 
So  the  Stepien  (1995)  results  should  also  be  roughly  applica¬ 
ble  to  our  system.  Given  the  main-sequence  lifetimes  of  our  pro¬ 
posed  components  are  «56  Gyr  and  the  age  of  the  Galactic  disk 
is  ~  12  Gyr  (Kaluzny  1990),  there  has  certainly  been  sufficient 
time  for  this  binary  to  reach  its  current  state  from  a  wider  initial 
separation.  Of  the  included  binaries  in  the  lists  of  near-contact 
binaries  (Shaw  1990,  Shaw  et  al.  1996),  very  few  have  spectral 
types  later  than  F  and  none  have  a  spectral  type  of  M,  making  this 
system  quite  rare. 

In  fact,  there  are  only  a  handful  of  stars  thought  to  be  of  sim¬ 
ilar  heritage.  YY  Gem  and  CM  Dra  (Strassmeier  et  al.  1988), 
CU  Cnc  (Ribas  2003),  BW3  V38  (Maceroni  &  Montalban 
2004),  and  RE  0618+75  (Jeffries  et  al.  1993)  are  potential  near¬ 
contact  siblings  along  with  their  contact  cousins,  for  example 
W  Cor  (Odell  1 996 )  and  VZ  Psc  ( Hrivnak  et  al.  1 995),  and  those 
found  listed  by  Hilditch  et  al.  (1988),  although  none  of  those 
contact  systems  have  a  dMe  component. 

It  is  likely  that  the  tidally  enforced  synchronous  rotation  of  the 
stars  in  this  system  helps  enhance  their  magnetic  activity.  Ro¬ 
tation,  especially  rapid  rotation,  is  a  prerequisite  for  the  creation 
of  strong  magnetic  phenomena  seen  in  stars  through  an  internal 
dynamo  (for  a  recent  review  see  Donati  2004).  FS  Aur-79  shows 
classic  signs  of  chromospheric  activity  through  the  emission  ob¬ 
served  in  Ho  and  H/3  from  our  spectroscopy  (Fig.  18)  and  star- 
spots  inferred  from  our  multiwavelength  light-curve  modeling.  It 
is  therefore  likely  that  this  binary  is  evolving  toward  contact  and 
has  been  for  a  long  time. 

The  temperatures  derived  by  light-curve  modeling  (7)  = 
4100  K,  73  =  3425  K)  are  supported  by  the  observed  spectra 
and  are  a  good  match  (Fig.  14)  for  a  dK7e  +  dM3e  near-contact 
system.  The  B  —  V  color  of  the  system  also  supports  these  spec¬ 
tral  types.  The  U  —  B  color  would  seem  to  be  inconsistent  with 
these  spectral  types,  indicating  a  hotter  primary  star.  However, 
because  of  the  presumed  starspots,  the  energy  output  in  the  UV 
is  likely  tainted  by  chromospheric  activity,  especially  in  the 


U  band  (Spruit  &  Weiss  1986;  Houdebine  et  al.  1996).  As  an 
independent  check  of  our  temperature  estimates,  the  depths 
of  TiO  absorption  bands  have  been  measured  as  described  by 
O’Neal  et  al.  (1998).  In  their  paper  they  obtained  empirical 
relationships  between  the  depth  of  TiO  absorption  bands  and 
temperature  (their  Fig.  1).  Measurements  of  the  TiO  bands  at 
7055,  7088,  and7127  Ain  our  high-resolution  spectra  are  con¬ 
sistent  with  our  temperatures  as  shown  in  Figure  19  for  TiO 
(7127  A). 

5.  CONCLUSIONS 

We  find  that  FS  Aur-79  is  a  near-contact  binary  system  at  a 
distance  of  ss  146  pc  with  dK7e  and  dM3e  components  separated 
by  1.62  Rq.  The  system  has  a  7-velocity  of +83  ±  5.5  km  s-1, 
with  the  primary  having  an  orbital  velocity  around  115  km  s_1 
and  the  secondary  predicted  to  have  an  orbital  velocity  near 
210  km  s_1.  The  primary  mass,  temperature,  and  radius  are 
«0.6  Mq,  4100  K,  and  0.67  R..  The  secondary  mass,  temper¬ 
ature,  and  radius  are  «0.3  M0, 3425  K,  and  0.48  Re.  This  system 
is  eclipsing  with  an  inclination  angle  of  83°,  with  unequal  eclipse 
depths  indicating  the  temperature  difference,  and  unequal  quad¬ 
rature  magnitudes  indicating  the  presence  of  starspots.  A  U  —  B 
excess  and  H  Balmer  lines  in  emission  are  further  evidence  of 
chromospheric  activity.  The  modeling  indicates  the  presence  of 
two  major  spot  regions  on  the  primary  and  one  on  the  secondary. 
It  is  likely  that  this  system  has  evolved  from  larger  separations 
through  angular  momentum  loss  via  stellar  winds  that  are  cou¬ 
pled  to  the  magnetic  fields  from  chromospheric  activity  and 
forced  to  corotate  out  to  large  radii,  thereby  extracting  angular 
momentum  from  the  system  and  tightening  the  orbit.  This  likely 
evolution  means  this  system  may  become  a  W  UMa  system  with 
unique  late-type  components,  which  is  rare  given  that  the  evolu¬ 
tion  time  is  so  long. 

Further  study  of  the  high-resolution  spectra  are  under  way, 
including  Doppler  tomography,  deconvolving  the  primary  and 
secondary  spectra,  and  detailing  of  the  components  of  emission 
features  as  a  function  of  phase.  This  should  provide  further  in¬ 
formation  concerning  the  starspots,  the  chromospheric  regions, 
and  refinement  of  the  stellar  masses. 
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